SECURITY  CLASSIFICATION  OF  THIS  PAGE 


m 


OCUMENTATION  PAGE 

1  lb.  RESTRICTIVE  MARKINGS 


Form  Approved 
OMB  No  07Q4-0188 


AD-A232  891 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S)  | 

N00014-89-J-1237 

6a  NAME  OF  PERFORMING  ORGANIZATION 

Colorado  State  University 

6b  OFFICE  SYMBOL 
(If  applicable) 

6c.  ADDRESS  (City.  State,  and  ZIP  Code) 

Department  of  Chemistry 
Fort  Collins,  CO  80523 

8a  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 

8b  OFFICE  SYMBOL 
(if  applicable) 

8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

000  North  Quincy  Street 
Arlington,  VA  22217-5000 

3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

ADoroved  for  public  release;  distribution 
uni imi ted . 


S  MONITORING  ORGANIZATION  REPORT  NUM3E»(Si 


1C  FILE  COPY 


7b  AOORESS  (City.  State,  and  ZIP  Code) 


N00014-39-J-1237 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 


PROJECT 

task 

NO 

NO 

1 1  TITLE  (Include  Security  Classification) 

Normal  Mode  Analysis  of  van  der  Waals  Vibrations 


12  PERSONAL  AUTHOR(S) 

in 


13a  TYPE  OF  REPORT  13b  TIME  COVERED 

Technical  Reoort  from _ To _ 


16  SUPPLEMENTARY  notation 


WORK  UNIT 
ACCESSION  NO 


COSATI  CODES 


GROUP 


SUB-GROUP 


19  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

SEE  ATTACHED  ABSTRACT  r  M  , 


DTIC 

I  ELECTR 
MAR  27 1801 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21  A8STRACT  SECURITY  CLASSIFICATION 

KXunclassifieo/unlimited  Q  same  as  rpt  □  DTic  users  Unci  a  s s  i  f  i  ed  _  _ 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  |22b  TELEPHONE  (Include  Area  Code)  22c  Office  SWSC. 


I/MIIH<41II 


DO  Form  1473,  JUN  86 


Previous  editions  are  obsolete. 

S/N  0102-LF-014-6603 


SECURITY  CLASS'F.CAr.Q.N  O-  PAG* _ 

91  3  20  037 


OFFICE  OF  NAVAL  RESEARCH 
Contract  N00014-89-J-1237 
TECHNICAL  REPORT  #72 


NORMAL  MODE  ANALYSIS 
OF  VAN  DER  WAALS  VIBRATIONS 


S.  Li,  E.  R.  Bernstein 

AcceE 

;icri  i  cr 

Submitted  to  the 

NTI 3  CRAv/.l 
DTiC  TA3 
Unaonoti.icoti 

JOURNAL  OF  CHEMICAL  PHYSICS 

i cation 

Chemistry  Department 

Colorado  State  University 

By 

Distribution  / 

' 

Av£j.:iabi;ii 

Fort  Collins,  CO 

March  7,  1991 

D 

A.  - 

Avu.i  • 
Sp* 

Reproduction  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government. 

This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited. 


ABSTRACT 


Two  alogrithms  are  presented  for  calculation  of  the  van  der  Waals  modes  of  weakly  bound 
clusters.  Both  methods  rely  on  the  harmonic  normal  mode  approximation  and  a  chosen 
intermolecular  potential.  These  calculational  techniques  differ  specifically  in  the  form  of  the  force 
field  employed:  one  method  uses  the  total  (both  intra-  and  inter-  molecular)  force  field  for  the 
cluster  and  the  other  uses  only  the  intermolecular  force  field.  Both  methods  require  rather  elaborate 
coordinate  transformation  and  their  first  and  second  partial  derivatives;  these  are  provided  in  detail. 
The  two  calculations  generate  eigenvalues  and  eigenvectors  that  are  in  complete  agreement  with  one 
another  for  a  given  potential.  The  methods  insure  that  the  van  der  Waals  modes  are  calculated  at 

the  proper  cluster  equilibrium  configuration  for  which  all  torques  and  forces  on  molecules  and/or 
atoms  are  effectively  zero.  Examples  are  presented  for  (H^O^  (n=  2,  ....  7),  benzene  clustered 

with  water,  methane,  and  ammonia,  and  a  number  of  different  intermolecular  potentials.  Some  of 
the  observed  heterogeneous  cluster  van  der  Waals  modes  are  reassigned  in  light  of  these  new 
results. 


I.  INTRODUCTION 

The  accurate  calculation  of  cluster  van  der  Waals  vibrational  modes  based  on  empirical 
potentials  is  of  central  importance  for  a  number  of  purposes:  1.  assignment  of  spectroscopically 
identified  cluster  vibrations;  2.  determination  and  validation  of  various  potential  energy  models  of 
cluster  interactions;  3.  determination  of  the  structure  of  both  large  and  small  clusters;  and  4. 
estimation  of  the  effect  of  the  mixing  of  molecular  internal  modes  with  cluster  intermolecular  van 
der  Waals  modes.  In  this  report  we  discuss  two  different  methods  for  the  proper  and  accurate 
calculation  of  van  der  Waals  modes  of  clusters:  one  method  calculates  van  der  Waals  modes  for 
fixed  molecular  geometries,  and  the  other  calculates  van  der  Waals  modes  following  a  distortion 
(relaxation)  of  molecular  geometry  brought  about  by  the  intermolecular  potential  energy.  We  refer 
to  these  algorithms  as  the  external  force  field  (EFF)  method  and  the  total  force  field  -  relaxed 
molecule  (TFF-RM)  method,  respectively. 

Systematic  study  of  intermolecular  motions  in  crystals  can  be  traced  back  to  the  early 
studies  of  Halford1,  Horning2,  and  Shimanouchi3.  Bernstein4  calculated  the  modes  of  a  benzene 
crystal  based  on  these  earlier  studies  and  Warshel  and  Lifson5  gave  general  formulas  for  crystal 
normal  mode  calculations  with  periodic  boundary  conditions.  Scheraga  and  co-workers  have 
presented  several  calculations6-7  of  van  der  Waals  modes  of  clusters.  More  recently  Menapace,  and 
Bernstein8  extended  this  approach  to  calculate  the  van  der  Waals  vibrational  modes  of 
heterogeneous  clusters  (i.e.,  CgHgCH^O)^  C6H6(Ar)j,  C6H6(CH4)j,  C6H6(NH3)j,  etc.).  Jortner 

et  al9-11  calculated  a  limited  subset  of  cluster  normal  modes  using  an  empirical  potential  and  an 

approximate  Schrddinger  equation.  Watts  et  al12'13  have  studied  the  structure  and  vibrations  of  the 

water  dimer  employing  a  number  of  different  approaches.  Janda  et  al14  have  studied  the  ethylene 

9-14 

dimer  and  ethylene  rare  gas  systems.  With  the  exception  of  the  quantum  mechanical  treatments 
the  above  studies  employ  the  GF  matrix  method  of  Wilson15  considering  both  inter-  and  intra¬ 
molecular  interactions  (force  fields)  simultaneously.  In  all  of  these  instances  the  original  molecular 
geometry  is  employed  in  the  calculation  -  any  molecular  distortion  due  to  the  intermolecular 


interaction  within  the  cluster  or  crystal  is  (implicitly)  ignored.  We  term  this  approach,  without 
molecular  structure  relaxation,  the  total  force  field-unrelaxed  molecule  (TFF-UM)  approach. 

Cluster  normal  mode  analyses  employing  the  GF  matrix  method  involve  (implicitly  or 
explicidy)  the  overall  cluster  translational  and  rotation  degrees  of  freedom.  The  TFF-UM  approach 
to  this  calculation  does  not  necessarily  set  these  modes  to  zero  energy,  as  of  course  a  correct 
calculation  must.  Indeed,  if  the  intermolecular  interactions  are  strong,  this  approach  gives 
significantly  large  cluster  overall  translational  and  rotational  energies  and  causes  extensive  mixing 
between  cluster  van  der  Waals,  translational  and  rotational  modes. 

Nonetheless,  the  GF  matrix  diagonalization  is  an  effective  method  for  determination  of  the 
harmonic  normal  modes  of  van  der  Waals  clusters;  how  to  generate  appropriate  cluster  F  (force 
field)  and  G  (reciprocal  mass)  matrices  is  the  concern  of  this  report.  We  discuss  two  algorithms 
for  the  generation  of  correct  G  and  F  matrices  for  cluster  normal  modes  of  vibration.  The  total 
force  field-relaxed  molecule  (TFF-RM)  method  involves  both  inter-  and  intra-molecular  force 
fields.  The  external  force  field  (EFF)  method  involves  rigid  molecules  and  the  intermolecular  force 
field  only.  The  results  of  both  calculations  are  essentially  identical  and  most  significantly  uniquely 
dependent  on  the  chosen  intermolecular  potential  only. 

In  the  ensuing  discussion  we  describe  the  potentials  and  concomitant  geometries  employed 
in  these  calculations.  General  descriptions  of  the  two  calculational  algorithms  (TFF-RM  and  EFF) 
are  then  presented.  In  the  third  section  the  details  and  manipulations  associated  with  both 
calculational  methods  are  given.  The  results  and  conclusions  are  presented  in  the  final  sections  of 
this  report. 

II.  GENERAL  DESCRIPTION  OF  THE  TECHNIQUES 
A .  Intermolecular  Potentials  and  Cluster  Geometries 

Cluster  intermolecular  potential  energy  is  calculated  from  several  empirical 
intermolecular  potential  forms.  Most  frequently  and  generally  the  atom-atom  Lennard- Jones  6-12- 
1  potential,  adapted  from  the  work  of  reference  16, 
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E  =  £  Ejj  =  S  [Ay  /  fy12  -  C(j  /  iy6  +Ciej/rij], 


(1) 


is  employed.  In  eq.  (1)  the  summation  is  over  all  atom  pairs  for  different  molecules,  parameters 
and  Cjj  are  given  in  reference  8,  r-  are  the  interatomic  distances  of  the  atoms  in  different  molecules, 
and  e-l  is  the  partial  change  on  atom  i.  The  can  be  obtained  from  the  work  of  Scheraga17'18  and 

also  from  MOPAC  519'20  calculation  using  an  appropriate  Hamiltonian. 

A  potential  employed  for  water  clusters  is  somewhat  different  and  is  adopted  from 
21 

the  work  of  Jorgensen: 


E  =  Xe  -  =  X  [d-./r.12  +  c./r..^  +  b-./r-  •■*]  +X'ee/r... 

i)  1  ij  ij  ij '  ij  y  ij  J  i  j  y 


(2) 


The  parameters  djj,  c^  and  bjj  are  given  in  reference  21.  The  summation  for  the 

Coulombic  term  in  the  potential  (denoted  by  X?  is  overall  the  charge  pairs  in  different  molecules, 

and  the  Lennard- Jones  summation  (denoted  by  X)  is  over  all  atom  pairs  for  atoms  in  o  fferent 

molecules.  Each  water  molecule  in  this  interaction  form  has  four  charges.  This  water  molecule  is 

thus  five  centered:  two  hydrogen  atoms  with  their  partial  charges  and  one  oxygen  atom  with  two 

partial  charges  located  at  approximate  lone  pair  positions.  This  potential  form  gives  the 

22 

experimentally  suggested  water  dimer  geometry. 

A  potential  form  proposed  for  benzene-water  clusters  is  an  atom-atom  Lennard- 
Jones  1-4-6-9-12  potential  adopted  from  the  work  of  Jonsson  et  al: 


E  =  SEij  =  I[As/,ij  +  Bij/ru4  +  Cij/r/  +  Du/^  +  Eu/ru12] 


(3) 


in  which  the  parameters  are  detailed  in  reference  23. 


Cluster  geometries  are  calculated  from  an  energy  minimization  routine  with  fixed 
(rigid)  molecule  structures:  in  the  routine,  the  molecules  comprising  the  cluster  are  randomly 
placed  in  space  and  the  molecules  are  translated  and  rotated  toward  the  cluster  structure  of  lower 
energy.  All  molecules  are  eventually  trapped  in  a  local  cluster  potential  energy  minimum.  An 
energy  minimum  geometry  is  defined  as  that  cluster  geometry  for  which  all  forces  and  torques  on 
a  molecule  are  less  than  ca.  10‘ 14  dynes  and  dyne-angstroms,  respectively,  and  all  molecular 
displacements  up  to  roughly  0.1  A  generate  higher  energy  cluster  structures. 

B .  Total  Force  Field  -  Relaxed  Molecule  Method 

This  approach  to  the  calculation  of  the  normal  modes  of  clusters  considers  both  the 
inter-and  intra-molecular  force  fields;  that  is,  the  total  cluster  force  field.  The  full  3n  x  3n  force 
field  matrix  F  in  Cartesian  coordinate  representation,  with  n  the  total  number  of  atoms  in  the 
cluster,  is  developed  and  employed. 

The  intramolecular  force  field  is  adopted  from  the  literature.  For  example,  the 
general  force  fields  for  ammonia,  water,  methane,  etc.  are  given  by  Herzberg24  and  the  general 
force  field  for  benzene  is  given  by  Whiffen.25  The  intermolecular  force  field  is  directly  calculated 
from  one  of  the  intermolecular  potential  energy  forms  (i.e.,  eq.  (1)  or  (3))  by  analytical  calculation 
of  second  partial  derivatives. 

In  order  to  obtain  meaningful  harmonic  force  constants  for  the  cluster  based  on  the 
analytical  second  derivatives  of  a  chosen  potential  form,  the  residual  force  (f)  on  all  the  atoms  must 
be  zero.  Since  the  energy  minimization  procedure  moves  rigid  molecules  as  a  whole,  the  force 
residue  on  each  atom  does  not  necessarily  get  minimized  as  the  force  residue  on  the  molecule  is 
minimized.  Thus  the  cluster  structure  is  no  longer  at  a  potential  minimum  when  the  internal  force 
field  is  subsequently  involved.  Force  residues  on  the  atoms  are  of  the  magnitude  10'5  to  10'6 
dynes,  compared  to  the  initial  force  and  torque  residues  on  each  molecule  of  10' 14  dynes  and 
dyne-angstroms,  respectively,  at  a  cluster  rigid  molecule  equilibrium  geometry.  Thus,  if  an 
accurate  reliable  harmonic  force  field  is  to  be  evaluated,  all  the  molecular  structures  in  the  cluster 
have  to  be  "relaxed"  in  accord  with  the  intra-  and  inter-molecular  potentials  employed  to  insure  that 


a  true  potential  minimum  for  the  atoms  in  the  cluster  has  been  achieved.  Higher  order  non¬ 
harmonic  corrections  to  these  intermolecular  force  fields  also  depend  essentially  on  the 
establishment  of  a  true  cluster  equilibrium  geometry  in  the  presence  of  all  forces,  torques  and 
interactions. 

In  the  calculations  discussed  below,  the  Newton-Raphson  method  is  used  to 
displace  each  atom  in  the  cluster  to  find  its  new  equilibrium  position.  Normal  mode  analysis  is 
applied  only  after  the  force  residue  on  each  atom  is  less  than  10*14  dynes,  van  der  Waals  mode 
energies  are  then  found  by  diagonalization  of  the  GF  matrix.  The  internal  molecular  modes  and  the 
six  cluster  translational  and  rotational  zero  energy  modes  are  then  identified  and  the  remaining 
modes  are  the  van  der  Waals  vibrational  modes  of  the  cluster.  Eigenvectors  are  calculated  to 
identify  and  confirm  the  character  of  all  the  generated  modes  of  the  cluster. 

The  quality  of  the  calculational  algorithm  can  in  part  be  judged  by  the  calculated 
"zero  energy"  translational  and  rotational  cluster  modes.  Of  course,  the  harmonic  assumption  is 
also  made  by  the  very  nature  of  the  calculation.  The  overall  quality  of  the  theoretical  results  for  van 
der  Waals  modes  of  clusters  (as  judged  by  comparison  with  experiments)  depends  on  this  above 
separation,  the  harmonic  approximation  and  on  the  quality  of  the  chosen  intermolecular  potential. 

C .  External  Force  Field  Method 

The  external  force  field  method  considers  each  (non-linear)  polyatomic  molecule  as 
a  rigid  body  with  only  six  degrees  of  freedom.  The  degrees  of  freedom  for  the  van  der  Waals 
clusters  are  reduced,  with  respect  to  those  needed  for  the  TFF-RM  method  first  discussed,  from 
3n-6  to  6m-6,  with  n  the  number  of  atoms  in  the  cluster  and  m  the  number  of  molecules  in  the 
cluster. 

Within  the  EFF  calculational  approach,  no  internal  force  field  is  involved.  The 
force  and  torque  residues  for  the  molecules  in  the  cluster  are  as  originally  calculated,  ca.  10' 14 
dynes  or  dyne-angstroms,  respectively.  The  mixing  between  the  van  der  Waals  modes  and  the 
cluster  translational  and  rotational  modes  caused  by  these  residual  first  derivative  terms  is 
negligible.  The  originally  obtained  cluster  geometries  can  be  used  directly  without  further 


modification.  The  6m  x  6m  GF  matrix  can  be  diagonalized  and  the  van  der  Waals  modes  of  the 
cluster  are  generated  along  with  the  six  zero  energy  modes  of  overall  cluster  translations  and 
rotations.  The  problem  posed  for  the  EFF  method  is  the  3n  to  6m  coordinate  transformation  and 
its  associated  first  and  second  order  derivatives:  these  are  discussed  in  detail  below. 

III.  ALGORITHMS  AND  COMPUTATIONAL  PROCEDURES 

In  this  section  the  formulas  used  in  both  the  TFF-RM  and  EFF  calculational  methods  are 
derived.  Several  coordinate  systems  are  required  for  these  calculations  and  they  will  be  briefly 
defined  first. 

A.  Coordinate  Systems  and  Coordinates 

We  will  use  and  refer  to  the  following  coordinate  systems:  1.  cluster  Cartesian 
coordinate  system;  2.  molecular  internal  coordinate  system,  involving  the  bond  lengths,  bending 
angles,  etc.  in  terms  of  which  the  molecular  internal  motions  are  described;  3.  a  molecule 
Cartesian  coordinate  system  with  its  origin  fixed  at  the  molecular  center  of  mass  and  its  axes 
parallel  to  those  of  the  cluster  Cartesian  coordinate  system;  and  4.  a  molecule  principal  axis  system 
(in  which  the  molecular  moment  of  inertia  tensor  is  diagonal).  Atoms  are  located  by  3  atomic 
coordinates  (x,  y,  z,)  and  molecules  are  located  by  six  coordinates  in  the  cluster  Cartesian  system. 

B .  Total  Force  Field-Relaxed  Molecule  Method 

A  function  g(x)  is  at  a  minimum  at  some  point  xQ  if  and  only  if  its  first  derivative 
dg(x)/8x  is  zero  and  its  second  partial  derivative  matrix  d2g(x)/dxidx.  (for  all  i,  j  =  1, ...,  n)  is 
positive  definite  at  this  point  xQ.  In  other  words,  a  true  cluster  potential  minimum,  with  zero  force 
residue  f(x)  =  dE/dx,  will  yield  a  semi-positive  definite  force  field  matrix  F  and  the  eigenvalues  of 
the  F  matrix  (or  the  GF  matrix)  will  thereby  be  all  positive  or  zero.  The  cluster  force  field  must  be 
calculated  at  the  equilibrium  position  of  each  particle  in  the  cluster  to  obtain  physically  meaningful 
(non-negative,  both  zero  and  positive)  normal  mode  energies. 

The  readily  available  and  adopted  forms  for  the  intermolecular  and  intramolecular 
interactions  and  force  fields  are  distinctly  different.  Intermolecular  interactions  are  represented  in 


the  form  of  a  potential  energy  function  dependent  on  the  atomic  position  (Cartesian)  coordinates. 
The  intermolecular  forces  and  the  intermolecular  force  fields  are  obtained  from  this  function  as  the 
sets  of  first  and  second  derivatives  (evaluated  at  the  equilibrium  configuration  chosen), 
respectively. 

On  the  other  hand,  the  intramolecular  force  fields  are  given  in  the  molecular  internal 
coordinate  system  in  terms  of  the  internal  displacement  coordinates.  The  molecular  internal 
potential  energy  and  restoring  forces  can  only  be  calculated  from  molecular  distortions  by 
assuming  the  molecular  internal  force  fields  to  be  harmonic.  Let  E,  f,  and  F  denote  cluster 
energy,  force  residue  vector,  and  force  field,  respectively.  Both  E  and  f  have  contributions  from 
the  intermolecular  (external)  interactions  and  the  intramolecular  (internal)  interactions,  such  that 


with 


and 


E,m  =  l  ASFiM(S)AS‘ 

Wx>  *  ASFim(S)«ASV3x). 


(4) 


In  eq.  (4)  the  subscripts  ext  and  int  denote  contributions  from  the  intermolecular 
interaction  and  the  molecular  distortions,  respectively,  AS  denotes  the  internal  displacement 
coordinate  vector,  x  or  S  in  parenthesis  indicates  the  coordinate  system  (cluster  Cartesian  or 
molecular  internal)  with  respect  to  which  the  matrix  or  function  is  presented,  the  x  subscript 
indicates  the  force  element  in  the  x-direction,  and  superscript  t  indicates  the  transpose  vector. 

The  force  field  F,  the  second  derivative  of  the  potential  E,  in  atomic  Cartesian 


coordinate  representation  can  be  written  as: 
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F(x)  =  Fext(x)  +  Fjnt(x) 

with 

Fint,xyW  =(3AS/3X)  Fim(S)  @ASV3y) 

+  AS  Fint(S)(92ASt/9xay).  (5) 

The  partial  derivatives  (dAS/dx)  and  (92AS/3x8y)  must  be  calculated.  Table  I  gives 
definitions  of  the  internal  coordinates  (bond  stretching,  angle  bending,  out  of  plane  wagging,  and 
torsion)  used  in  this  work.  Others  can  be  generated  as  required.  Most  of  these  definitions  appear 
in  reference  15,  but  a  few  of  the  coordinates  are  redefined  in  order  to  remove  potential 
singularities. 

The  first  partial  derivatives  of  molecular  internal  coordinates  with  respect  to  atomic 
Cartesian  coordinates  (expressed  in  either  the  cluster  or  molecule  Cartesian  coordinate  system  -  the 
elements  of  the  "B-matrix")  are  deduced  and  tabulated  in  reference  15.  The  analytical  second 
partial  derivatives  of  the  molecular  internal  coordinates  with  respect  to  atomic  Cartesian 
coordinates  have  not  been  previously  presented.  While  these  derivatives  are  straightforward  to 
obtain,  their  presentation  in  vector  form  is  useful  for  computations:  the  full  deduction  of  the 
requisite  first  and  second  partial  derivatives  in  vector  form  is  presented  in  Table  II. 

The  calculation  procedure  is  as  follows:  1.  cluster  energy  E,  force  residue  f,  and 
force  field  matrix  F  are  calculated  in  atomic  Cartesian  coordinate  representation  for  a  given 
geometry;  2.  if  the  magnitude  of  f  is  greater  than  10'14  dynes,  a  Newton-Raphson  method  is 
employed  to  find  atomic  displacements  toward  a  lower  cluster  potential  energy;  3.  the  total 
molecular  distortion  is  updated  in  molecular  internal  coordinates  for  each  molecule  of  the  cluster; 
4.  E,  f,  F  are  recalculated  as  above;  5.  termination  occurs  when  the  magnitude  of  f  is 
appropriately  small  (ca.  10' 14  dynes)  and  the  current  F  matrix  is  then  used  to  construct  the  GF 
matrix  product  for  normal  mode  analysis. 
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Note  that  the  molecular  distortion  in  internal  coordinates  AS  must  be  explicitly 
calculated  and  saved  in  order  to  calculate  the  molecular  restoring  force.  AS  cannot  be 
approximated  as  (3AS/3x)*Ax  as  the  atomic  displacements  become  large  (ca.  0.1  A). 

An  eigenvalue  shift  routine  is  also  employed  to  ensure  rapid  approach  to  a  true 
potential  energy  minimum  position.  Care  must  be  taken  to  reduce  computer  rounding  errors  in 
these  calculations  because  the  requirements  on  zero  cluster  translational  and  rotational  residual 
forces  are  quite  stringent.  These  force  residues  are  made  to  go  to  zero. 

C .  External  Force  Field  Method 

In  the  external  force  field  method,  only  the  intermolecular  interaction  energy  is 
considered;  that  is,  E  =  E  .  The  cluster  intermolecular  potential  energy  is  given  originally  as  a 

function  of  the  cluster  Cartesian  coordinates  x  of  the  atoms  and  any  charges  not  located  at  atomic 

positions.  The  intermolecular  interaction  energy  must  thus  be  transformed  into  a  function  of 
molecular  coordinates  Qj  for  the  jth  molecule  of  the  cluster  -  the  index  j  runs  over  all  m  molecules 

of  the  cluster.  In  order  to  get  both  F  and  G  matrices  in  this  molecular  displacement  coordinate 
representation, 

E  =  EextW  =  EextWQ))’  <6> 

the  appropriate  coordinate  transformation  must  be  obtained.  The  vector  Q  is  of  the  form  Q  =  (Qj, 
Q2.  with  Qj  given  below. 

1.  The  G  matrix15 

For  each  molecule,  the  molecular  displacement  coordinates  involve  three 
rotations  (<J>a,  <j>p,  4>^)  along  the  three  principal  axes  of  the  molecules  and  three  translations  (X,  Y, 

Z)  along  the  cluster  Cartesian  coordinate  system  axes.  One  can  then  write: 
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Qj  =  (Xj.Yj,Zj.VVV 

dXj  =  Z  mjdXj  /  Znv 

%  =  2m,  (7) 

in  which  nv  and  are  the  mass  and  x  coordinate  of  atom  i  (in  a  particular  molecule  j),  Laa  is  the 

molecule  moment  of  inertia  along  the  a-axis,  and  (Xj,  are  the  atomic  Cartesian  coordinates  of 

the  i 111  atom  in  the  jth  molecule  principal  coordinate  system.  Cyclic  permutation  of  the  variables 
produces  the  other  required  relations.  Based  on  these  relations,  the  molecular  B-matrix  (dQJdx) 

can  now  be  formed.  The  Gj  matrix  is  thereby  given  as,15 

Gj =  Bj  Mj"1  Bj1,  (8) 

in  which  NT  is  the  (diagonal)  molecular  atomic  mass  matrix.  Eq.  (8)  yields  the  correct  diagonal 
molecular  G  matrix,  with  the  inverse  of  the  molecular  mass  as  the  entry  for  the  molecular 
translational  motions,  and  the  inverse  of  the  corresponding  moment  of  inertia  for  the  molecular 
rotational  motions. 

Since  the  translational  and  rotational  degrees  of  freedom  of  one  molecule  are 
independent  of  those  of  another,  the  molecular  displacement  coordinates  of  one  molecule  are 
orthogonal  to  those  of  another.  The  B  matrix  of  the  whole  cluster  is  thus  block  diagonal  and  the 
G  matrix  of  the  cluster  is  diagonal. 

2.  The  F  matrix15 

To  find  the  force  field  matrix  F  in  the  same  coordinate  system  we  return  to 

eq.  (6)  to  find 


9E/8Q  =  0E/8x)*0x/8Q) 


(9) 


and  one  element  of  F(Q)  for  a  given  molecule  is, 
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32E/3QvaQ4  =  0x/aQv).02E/ax3y).0y/aQ4)  +  0E/3x)*O2x/3Qv8Cy  (10) 

in  which  Q,  Q  or  Q  indicates  any  of  the  molecular  coordinates  X;,  Y-,  Z.,  <J>.  ,  01R,  <J>.  for 
v  h-  j  j  j  jp  jy 

all  j. 

Since  8E/9x  is  the  force  f  on  atoms  and  not  necessarily  zero  in  the  EFF 
approach,  both  (dx/dQ)  and  (92x/8Qv8Q^)  are  required  to  calculate  F(Q).  Note  that  the  (dx/dQ) 

cannot  be  obtained  from  inversion  of  the  B  matrix  as  B  is  singular.  The  analytical  determination 

of  these  first  and  second  derivatives  is  discussed  in  the  ensuing  presentation. 

Atomic  displacements  dx  in  one  molecule  are  independent  of  molecular 
displacements  dQ  of  another  molecule  and  the  matrices  (dx/dQ)  and  (d2x/dQvdQ^)  are  thereby 

block  diagonal  in  terms  of  molecules  in  the  cluster.  We  thus  drop  any  notation  and  indices 
referring  to  a  particular  molecule. 

For  a  rigid  molecule,  an  atomic  transitional  displacement  is  just  the 
molecular  translational  displacement  and  thus, 

dx-JdX  =  1 ,  dxi  Id  Y  =  0  and  d\/dX  9Q^  =  0  (11) 

in  which  Xj  are  atomic  coordinates  of  the  ith  atom,  X  and  Y  are  the  molecular  x-  and  y- 
coordinates,  respectively,  and  Q  can  be  any  molecular  translational  or  rotational  coordinate. 

r* 

Other  derivatives  are  found  by  cyclic  permutation  of  the  coordinates  in  eq.  (1 1). 

The  derivatives  of  eq.  (9)  and  (10)  must  be  found  only  for  the  molecular 
rotational  coordinates  (<J>a,  <J>p,  <f>^).  The  derivatives  can  be  explicitly  generated  by  first 

considering  the  general  rotation  matrix  R(0)  for  a  rotation  about  an  arbitrary  direction  in  space. 
Let  (a},  ^  a3)  be  a  unit  vector  along  that  direction  and  let  0  be  the  rotation  angle  about  this 

direction.  The  rotation  matrix  R(0)  is  then, 26 
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R(0)  = 


^cos0  +  (1  -  cos0) 

a^  (1  -  cos0)  +  a^in© 
^a^  (1  -  cos0)  -  a^in© 


ata2(l  -  cos0)  -  a3sin0 
cos0  +  a^  (1  -  cos0) 
a^Q  -  cos0)  +  ajSin© 


ata3(l  -  cos0)  +  a2sin0  ^ 
a2a3(l  -  cos0)  -  ajSin© 
cos0  +  a^(l  -  cos0) 


(12) 

This  rotation  matrix  is  derived  as  a  product  of  three  successive  rotations:  a 
rotation  of  the  z  axis  to  the  direction  (a^  a2>  a3);  a  rotation  of  the  coordinate  system  around  the 
new  z  axis  (a^  ^  a3)  by  an  angle  -0;  and  a  rotation  of  the  z  axis  back  to  its  original  direction. 

Based  on  this  general  rotation  matrix  form,  an  ordered  set  of  rotation 
matrices  Ra  (<]>a),  Rp(<J)p),  RyC^y)  f°r  each  atom  (i)  in  the  molecule  can  be  defined  as, 

r  =  Ra(<t,a)Rp(<l,p)Ry(4,y)ro  • 

In  eq.  (13)  Ra,  Rp,  and  R^  are  the  matrix  representations  of  rotations  along 
the  molecular  principal  axes  a,  (5,  and  y,  respectively,  rQ  are  the  unrotated  and  r  the  rotated 

atomic  coordinate  vectors  in  the  molecular  Cartesian  coordinate  system.  With  this  definition,  we 
have  for  each  atom  in  a  molecule, 

3rV9<t>a  =  dRa(<J>a)/d(t)aRp((t)p)Ry(4)y)r^ 

d2rl/d$l  =  d2Ra(4>a)/d4)2aRp(4)p)R7(4)y)r‘ 
a2rt/a0cta<J)p  =  dRa^yd^dRp^Vd^R/^rl, . 

When  evaluated  at  4>a  =  <t>p  =  $y  =  0.  these  derivatives  are, 


(14) 
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ar7a<j>a = 


a  V/9<t£  = 


0 

"ai3 

ai2  ^ 

ai3 

0 

-ail 

To 

"ai2 

all 

0  J 

aiiai2 

aiiai3 

aiiai2 

4-i 

ai2ai3 

^1% 

a!2a!3 

4-i 

J 


(15a) 


(15b) 


a2rt/a<t)0La<t)p  = 


■iAi 


ai2a21 

ai2a22 

ai2a23 


ai3a21 

ai3323 


=  an  (aa  •  f0) 


The  other  derivative  forms  not  explicitly  stated  in  eq.  (15)  can  be  obtained 
by  cyclic  permutation  of  the  indices.  In  these  equations  aa  =  (an,  a12,  a13),  ap  =  (a21,  a^,  a^) 

and  a^=  (a31,  a32,  a33).  The  ay  are  the  three  unit  vectors  along  the  three  molecule  principal  axes 

represented  in  the  original  cluster  Cartesian  coordinate  system.  Their  components  are  simply  the 
direction  cosines  of  the  molecular  principal  axes  in  the  cluster  Cartesian  coordinate  system  if  the 
rotation  is  a  proper  rotation.  The  rotation  matrix  formed  by  these  vectors, 


P  =  (aoi,  a^,  a^) 


<*2! 

ai2 

*32 

*33 

(16) 


can  be  gotten  from  the  orthogonal  matrix  which  diagonalizes  the  molecular  moment  of  inertia 


tensor. 
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One  can  readily  demonstrate  that  the  rotation  order  (order  of  the  matrices  in 
eq.  (13))  is  of  no  consequence  for  the  final  F  (Q)  provided  that:  the  same  rotation  order  is 
maintained  for  all  atoms  in  a  molecule;  the  molecule  has  no  net  torque;  and  the  rotation  is  an 
infinitesimal  one.  Moreover,  the  rotation  matrix  P  in  eq.  (16)  may  be  either  proper  or  improper 
(i.e.,  IPI  =  ±  1).  Of  course,  all  the  rotation  matrices  for  all  the  molecules  in  the  cluster  must  be 
simultaneously  proper  or  improper. 

3.  Eigenvectors 

To  visualize  the  resulting  eigenvectors  of  the  GF  matrix  in  the  EFF  method, 
transformation  of  these  eigenvectors  into  the  cluster  Cartesian  coordinate  system  is  quite  helpful. 
We  discuss  this  procedure  below. 

Three  successive  small  rotations  (0a,  0^,  0^)  of  a  molecule  around  its  three 
molecular  principal  axes  are  equivalent  to,a  single  rotation  0  around  the  axis  pointing  in  the 

0p> 

through  (0a,  0^,  0^)  in  the  molecular  principal  coordinate  system  is  equivalent  to  rotation  of  the 
molecule  about  (0a,  0p,  0^)  P  in  the  molecular  coordinate  system  with  axis  parallel  to  those  of  the 

original  cluster  Cartesian  coordinate  system. 

For  three  infinitesimal  rotations  0a,  0p,  0Y  around  the  three  molecular 

principal  axes,  eq.  (15a)  gives. 


V 


Rotation  of  a  molecule 


direction  (0a,  0p,  0^)  for  which  0  is  the  norm  of  the  vector  (0a, 


rl  =  (r0  +  dr)1  = 


1 

0'T  1 
(-fu  f'o 


f’p ) 

•♦'a 


1 


J 


(17) 


in  which  (0’a,  0'p,  0'y)  =  (0a,  0p,  0y)  P. 


Let  the  eigenvector  coefficients  for  a  molecule  j  be  (XL,  YV,  Z L,  <{)ja,  <f)jp,  0^.  The  equivalent  but 
unnomalized  eigenvector  coefficients  for  each  atom  i  in  molecule  j  are, 
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(  0 


V 


-V 


V 

-V 


(x:\ 

j 

r!  + 

y; 

j 

1 

/ 

X' 

^  J  / 

(18) 


in  which  again  (<j>'ja,  (j)^,  <>' V)  =  (<f>ja,  4>jp,  <J>^)  Pj,  and  Tj  is  the  atomic  coordinates  of  atom  i  in  the 

molecular  coordinate  system  with  axes  parallel  to  those  of  the  original  cluster  Cartesian  coordinate 
system. 

4 .  Summary  of  EFF  Method 

The  EFF  computational  procedure  is  summarized  as  follows:  1.  calculate 
the  intermolecular  force  vector  f  and  force  field  F  in  atomic  coordinates  based  on  the  cluster 
geometry  of  the  local  energy  minimum  with  rigid  molecular  structure;  2.  find  the  molecular 
moment  of  inertia  tensor  for  every  molecule  in  the  cluster  and  diagonalize  it;  3.  the  molecular  mass 
and  principal  moments  of  inertia  are  saved  as  the  inverse  entry  of  the  usual  G  matrix;  4.  the 
orthogonal  matrices  which  do  this  diagonalization  are  determined  to  be  proper  or  improper  and  are 
saved  as  the  rotation  matrices  P  in  eq.  (16);  5.  the  matrices  of  the  first  partial  derivatives  9x/9Q 
and  second  partial  derivatives  92x/9Qy9Q^  are  calculated  based  on  eq.  (1 1)  and  (15);  6.  the  GF 

matrix  is  formed  in  the  molecular  coordinate  representation  and  diagonalized;  and  7.  the 
eigenvectors  are  transformed  back  to  atomic  coordinates  for  visualization  according  to  eq.  (18). 

IV.  RESULTS  AND  DISCUSSION 
A.  Homogeneous  Clusters 

Table  III  contains  the  results  of  calculations  of  the  van  der  Waals  modes  of  the 
water  dimer  by  several  different  techniques:  TFF-UM,  TFF-RM,  EFF,  and  MOPAC  5  algorithms. 
Different  potential  forms  are  employed  as  well  for  the  two  methods  advocated  in  this  work. 
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Unfortunately,  no  experimental  results  are  available  for  comparison.  The  cluster  geometries 
associated  with  these  potentials  are  given  in  Figure  1.  The  "relaxed"  geometry  obtained  for  this 
dimer  (TFF-RM)  is  not  distinguishable  from  that  of  the  original  one  presented.  The  results  for  the 
two  methods  discussed  in  detail  (TFF-RM  and  EFF)  are  virtually  indistinguishable  -  this 
demonstrates  the  detailed  correctness  of  both  calculational  methods.  The  residual  non-zero 
rotational  and  translational  energies  are  due  almost  entirely  to  round-off  errors  in  the  computer 
calculation. 

Since  no  five-center  water  molecule  internal  force  field  is  available,  mixing  of  the 
eq.  (2)  intermolecular  force  field  and  any  three-center  intramolecular  force  field  in  any  TFF 
approach  causes  non-zero  cluster  translational  energies.  This  potential  then  can  only  be  employed 
within  the  EFF  method. 

The  TFF-UM  approach  with  any  intra-  and  inter-  molecular  potential,  generates 
large,  residual  forces  on  the  atoms  and  thereby  large  non-zero  (negative)  eigenvalues  for  the  overall 
translational  and  rotational  modes  of  the  clusters.  This  is  seen  in  Table  III. 

One  might  anticipate  that  by  increasing  the  internal  force  field  and  thereby  making 
the  molecule  more  rigid,  one  could  force  these  negative  eigenvalues  to  zero.  The  results  presented 
in  Table  IV  do  not  support  this  position.  Two  reasons  can  be  cited  for  the  failure  of  this  approach: 
1 .  the  force  residue,  from  which  these  negative  eigenvalues  derive,  is  neither  reduced  nor  balanced 
by  increasing  the  internal  force  field;  and  2.  the  disparity  between  the  sizes  of  the  two  force  fields 
eventually  causes  the  matrix  diagonalization  routine  to  fail. 

The  first  column  of  Table  IV  gives  the  calculated  energies  for  the  water  dimer 
geometry  given  in  Figure  lc.  This  structure  is  at  a  saddle  point  on  the  intermolecular  potential 
surface  described  in  eq.  (1).  Cluster  geometries  at  true  energy  minima  (f  =  0  and  IFI  semi-positive 
definite)  must  have  six  zero  energy  translational  and  rotational  modes,  and  no  negative  vibrational 
energies.  Cluster  geometries  at  saddle  points  may  also  have  six  zero  energy  modes  but 
additionally  have  negative  energy  vibrational  modes  due  to  the  negative  curvature  of  the  potential 
surface  at  the  saddle  point  in  some  direction.  Non-zero  residual  forces  will  cause  non-zero 
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translational  and  rotational  energies.  Thus  a  van  der  Waals  mode  analysis  can  be  used  to  check 
for  true  minimum  energy  or  saddle  point  cluster  structures.  These  results  (Table  III  and  IV) 
suggest  that  the  geometries  observed  from  the  energy  minimization  analysis  are  indeed  at  local 
energy  minima. 

Larger  water  clusters  [(I^O^,  n  =  3, ...,  see  Figure  2]  represent  a  very  difficult 

calculation  for  the  TFF-UM  approach  due  to  the  strong  intermolecular  interactions.  Thus,  mixing 
of  the  van  der  Waals  modes  with  both  the  internal  modes  of  water  and  the  overall  rotational  and 
translational  motions  of  the  cluster  (f  >  0)  can  be  quite  extensive.  Table  V,  along  with  Tables  III 
and  IV,  provides  an  example  of  this  problem.  Clearly  the  normal  modes  of  even  the  simplest 
clusters  can  be  reliably  calculated  only  after  the  residual  forces  in  the  cluster  are  zeroed  (by  TFF- 
RM,  EFF,  or  some  other  technique). 

The  difference  between  mode  energies  calculated  by  the  TFF-RM  and  EFF  methods 
is  greatest  for  the  high  energy  modes  due  to  the  mixing  of  the  van  der  Waals  modes  with  the 
internal  modes:  the  TFF-RM  high  energy  modes  are  lower  in  energy  than  the  comparable  EFF 
modes  as  is  to  be  expected.  The  two  methods  yield  nearly  identical  results  for  the  low  energy 
modes.  The  requisite  relaxed  geometry  calculated  for  the  TFF-RM  method  is  different  from  that 
of  the  original  energy  minimization  result  by  typically  less  than  0.1  A  (see  Figure  2). 

B .  Heterogeneous  Clusters 

van  der  Waals  energies  calculated  by  the  TFF-UM,  TFF-RM,  and  EFF  methods  for 
various  different  potentials  (eq.  (1)  and  eq.  (3))  are  presented  in  Table  VI  for  C6H6(H20)j, 
(NH3)1  and  (CH4)j.  These  results  are  compared  to  the  earlier  calculations  of  reference  8.  The 

maximum  cluster  translational  and  rotational  energies  for  each  calculation  are  given  in  parenthesis. 
All  the  corresponding  cluster  geometries,  based  of  course  on  the  various  potentials,  are  depicted  in 
Figure  3. 

The  calculational  results  for  heterogeneous  clusters  presented  in  Table  VI  are  nearly 
identical  for  both  the  TFF-RM  and  EFF  methods.  Apparently  little  internal  mode-external  mode 
coupling  takes  place  in  these  clusters.  Clearly  cluster  geometry  (i.e.,  C6H6(NH3)j)  and  cluster 
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intermolecular  potential  are  the  major  determining  factors  for  normal  mode  energies  in  these 
clusters  if  the  calculations  are  performed  at  an  appropriate  equilibrium  position. 

Based  on  these  new  calculations  we  can  attempt  to  reassign  the  observed  van  der 
Waals  modes  for  the  clusters  presented  in  Table  VI  and  Figure  3.  The  right-hand  column  contains 
the  experimentally  observed  van  der  Waals  mode  energies.  Only  the  fundamentals  that  seem 
reasonably  intense  and  well  characterized  are  included  in  this  listing.  If  the  molecular  structures 
are  not  properly  relaxed  (TFF-UM)  to  drive  f  to  zero,  the  low  frequency  van  der  Waals  modes 
mix  with  the  cluster  translational  and  rotational  motions.  The  extent  of  mixing  not  only  depends 
on  the  magnitude  of  the  force  residues,  but  also  on  the  nature  of  the  internal  force  field.  This  latter 
point  can  be  seen  in  the  first  two  columns  of  Table  VI  for  the  TFF-UM  calculation:  the 
comparison  is  for  the  data  of  reference  8  and  the  calculations  presented  herein  using  a  more 
complete  and  accurate  force  field  for  both  molecules  of  the  heterogeneous  cluster. 

The  eigenvectors  calculated  by  both  the  TFF-RM  and  EFF  methods  are  identical  for 
a  given  potential  function;  the  form  of  the  eigenvectors  and  the  mode  energies  are,  of  course, 
highly  dependent  on  the  chosen  intermolecular  potential.  This  can  be  seen  in  Figure  4  quite 
clearly. 

C .  TFF-RM  or  EFF? 

While  the  EFF  method  of  calculating  cluster  normal  modes  of  vibration  does  not 
yield  the  internal-external  mode  mixing  and  concomitant  mode  splittings  and  shifts,  it  does 
provide  a  simple  effective  method  that  yields  quite  reasonable  results,  especially  if  the  external  and 
internal  modes  are  well  separated  in  energy.  For  large  clusters  we  recommend  the  EFF  method 
(with  the  above  caveats  )  for  the  following  reasons:  1.  only  the  intermolecular  potentials  are 
required;  2.  reminimization  of  the  cluster  energy  is  not  required;  3.  since  each  molecule  is 
described  by  only  six  coordinates  much  less  computer  memory  is  needed  for  this  method 
compared  to  the  TFF-RM  method;  and  4.  the  TFF-RM  method  can  be  applied  only  to  atom-atom 
interaction  potentials  with  point  charges  located  at  atomic  mass  positions. 


21 


V.  CONCLUSIONS 

Two  methods  for  the  calculation  of  van  der  Waals  modes  of  molecular  clusters  are 
presented.  One  uses  both  internal  and  external  force  fields  to  calculate  all  the  vibrational 
eigenvectors  and  eigenvalues  for  a  given  cluster.  The  other  employs  the  intermolecular  force  field 
only  and  thus  obtains  only  van  der  Waals  mode  eigenvectors  and  eigenvalues.  Both  techniques 
ensure  that  modes  are  calculated  for  equilibrium  structures  of  the  cluster  only.  Both  approaches 
give  excellent  separation  between  van  der  Waals  modes  and  overall  cluster  (zero  energy) 
translational  and  rotational  degrees  of  freedom  and  give  no  negative  eigenvalues. 

The  TFF-RM  method,  employing  both  internal  molecular  and  intermolecular  force  fields, 
shows  small  mixing  between  the  lower  energy  molecule  modes  and  the  high  energy  cluster  modes. 
The  two  methods,  TFF-RM  and  EFF,  yield  virtually  identical  results  for  van  der  Waals  modes 
below  150  cm'1,  for  the  clusters  and  potential  energy  forms  considered  herein.  In  general,  the 
calculation  seems  to  be  in  good  agreement  with  the  sparse  experimental  data.  The  TFF-RM 
method  also  generates  "exchange  or  exciton"  type  interactions  between  the  internal  molecule  modes 
of  the  clusters. 
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Table  I.  Definitions  of  Internal  Coordinates 

Internal  Mode  Variable,  definition,  and  Formula  Expression 

Bond  Stretching  rjj  (distance  between  atom  i  and  j) 

rij  =  Injl  =  In  -  rjl 

Angle  Bending  I  0}jk  (angle  between  bond  r jj  and  rkj) 

Cos  (9 ijk)  =  **ij»r kj/rijrk j 


Angle  Bending  II  0h  (angle  created  by  bond  rjj  bending  with  respect 

to  the  bisector  of  the  k-j-1  angle)  0h  =  (0ijk  -  © iji)/2. 

Out  of  plane  Yijki  (angle  between  rjk  and  the  plane  determined  by 
Wagging  j,  k  and  1  atoms) 

Sin  (yijki )=  r,k*(rjk  x  rki)/riklrjk  x  rkjl 


Torsion3  d>ijkl  (angle  between  bonds  rjj  and  rki,  viewed  along 

rjk  bond) 

Sin(Oijki)  =  ((rjj  x  rjk)  rk|)  rjk  /  Irjj  x  rjkl  lrjk  x  rk|l. 


aThis  definition  of  d>jjki  is  valid  only  about  O  ~  0  and  ~  jt. 


At  <h  ~  —  and  ~ 
2 


3k 

2 


this  definition  yields  a  singularity  for  the  necessary  derivatives  and  a  new 
relation  is  required. 


Table  II.  Derivatives  of  Vector  Quantities 
a,  b,  c,  d  are  vectors,  and  a,  b,  c,  d  are  their  lengths. 


Derivatives  of  u  =  a*b/(aM 

3u/3x  =  [3(a*b)/3x]/(ab)  -  (3a/3x)u/a  -  (3b/3x)u/b 

32u/3x3y  =  [a2(a»b)/3x3y]/(ab) 

-  [(3a/3y)/a  +  (3b/ay)/b](au/3x)  -  [(aa/ax)/a  +  (ab/ax)/b](au/ay ) 

-  [(aa/ax)(ab/ay)  +  (3a/3y)(ab/3x)]  u/(ab) 

-  (32a/3x3y)  u/a  -  (32b/3x3y)  u/b 
Derivatives  of  v  =  (a»d)c/(ab) 

3v/3x  =  [3(a*d)/3x]  c/(ab)  +  (a*d)(3c/3x)/(ab) 

-  (3a/3x)v/a  -  (3b/3x)v/b 

32v/3x3y  =  [32(a«d)/3x3y]  c/(ab)  +  (a«d)(32c/ax3y)/(ab) 

+  {[a(a«d)/ax](3c/ay)  +  [3(a«d)/3y](3c/3x))/(ab) 

-  [(3a/3y)/a  +  (3b/ay)/b](3v/ax)  -  [(3a/3x)/a  +  (3b/ax)/b](3v/ay ) 

-  [(3a/3x)(3b/3y)  +  (3a/3y)(ab/ax)]  v/(ab) 

-  (32a/3x3y)  v/a  -  (32b/3x3y)  v/b 

Other  Necessary  Derivatives 

3a/3x  =  a»(3a/3x)/a 

32a/3x3y  =  [  a*(32a/3x3y)  +  (3a/ax)*(3a/ay)  -  (3a/3x)(3a/3y)]/a 

3(a»b)/3x  =  (3a/3x)»b  +  a»(3b/3x) 

3(a  x  b)/3x  =  (3a/3x)  x  b  +  a  x  (3b/3x) 

When  both  a  and  b  are  ry : 

32(a»b)/3x3y  =  (aa/3xHab/3y)  +  (3a/ay)*(3b/ax) 

32(a  x  b)/3x3y  =  (3a/3x)  x  (3b/3y)  +  (3a/3y)  x  (3b/3x) 


Table  III.  vdW  normal  mode  comparisons  (in  cm-1)  for  (H20)2. 
The  entry  in  parenthesis  in  each  column  represents  the  largest 
non-zero  energy  for  the  cluster  translational  and  rotational 
motions. 


Eq.  1 

Eq.  1 

Eq.  2 

MOPAC  5 

Eq.  1 

TFF-RM 

EFF 

EFF 

TFF-UM 

441.4 

459.0 

723.1 

483.5 

437.3 

289.5 

291.0 

716.6 

434.1 

238.5 

153.5 

155.8 

307.5 

379.6 

108.5 

1 18.5 

119.4 

301.1 

270.6 

57.5 

95.9 

96.9 

263.4 

190.7 

0.2 

26.9 

26.9 

99.6 

100.2 

0.0 

(0.7) 

(0.1) 

(0.2) 

(121.0) 

0.0 

-17.5 

-33.0 

-54.5 

-69.7 

-85.6 


Table  IV.  The  effect  of  a  saddle  point  structure  and  an  increased 
internal  force  field  on  cluster  normal  modes  (in  cm-1):  (H20)2- 
Example  for  TFF-UM  employing  Eq.  1. 


Saddle  Point 

Field  x  10° 

Field  x  102 

Field  x  104 

Field  x 

438.5 

437.3 

456.3 

458.1 

1126.7 

245.8 

238.4 

238.5 

249.2 

702.3 

103.2 

108.5 

109.1 

133.7 

625.6 

97.7 

57.5 

57.5 

57.5 

395.5 

90.6 

0.2 

0.0 

43.5 

61.9 

0.1 

0.0 

0.0 

0.0 

57.3 

0.1 

0.0 

-1.9 

0.0 

0.0 

0.0 

-17.5 

-17.3 

-14.7 

0.0 

0.0 

-33.0 

-33.0 

-23.8 

-14.8 

0.0 

-54.5 

-54.4 

-46.7 

-24.5 

0.0 

-69.7 

-69.7 

-56.9 

-52.0 

-409.2 

-85.6 

-85.5 

-79.0 

-66.2 

TABLE  V.  van  der  Waals  Mode  Energies  (cm'1)  for  a  Few  Water  Clusters 
employing  the  Intermolecular  Potential  given  in  Eq.  1 
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TABLE  V.  (cont’d) 
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Table  VI.  Normal  mode  (in  cm-1)  comparisons  for  inhomogeneous 
clusters.  The  structures  are  depicted  in  Figure  3.  The  entry  in 
each  column  in  parenthesis  represents  the  largest  value  for  the 
set  of  modes  indicated. 


C6H6(H20)1 


Ref.  8 

Eq.  1 

Eq.  1 

Eq.  1 

Eq.  3 

Expt. 

Mode  Identity 

TFF-UM 

TFF-RM 

EFF 

EFF 

(Approx.) 

156 

167.5 

162.7 

163.0 

287.0 

Tx 

50 

163.4 

161.4 

162.5 

125.0 

Ty 

159 

51.3 

50.2 

50.2 

82.0 

Stretch 

18 

30.0 

29.5 

29.4 

39.0 

30. 3a 

Bendl 

14 

21.1 

24.1 

24.2 

38.8 

25. 5a 

Bend2 

40 

40.0 

4.4 

4.3 

5.1 

5.2b 

Tz 

(?) 

(6.31) 

(0.49) 

(0.04) 

(0.05) 

(Rot.) 

(?) 

(0.44) 

(0.00) 

(0.00) 

(0.00) 

(Tr.) 

C6H6(NH3)i-Geometry  A 


Ref.  8 

Eq.  1 

Eq.  1 

Eq.  1 

Mode  Identity 

TFF-UM 

TFF-RM 

EFF 

Expt. 

(Approx.) 

152 

156.1 

152.0 

153.1 

Tx 

152 

156.1 

152.0 

153.1 

Ty 

97 

97.7 

97.8 

98.4 

Stretch 

19 

30.7 

29.5 

29.7 

Bendl 

19 

30.7 

29.5 

29.7 

Bend2 

44 

43.2 

16.0 

16.3 

15. 0b 

Tz 

(?) 

(7.96) 

(0.65) 

(0.06) 

(Rot.) 

(?) 

(0.41) 

(0.00) 

(0.00) 

(Tr.) 

Table  VI  (continued) 


C6H6(NH3)i-Geometry  B 


Ref.  8 

Eq.  1 

Eq.  1 

Eq.  1 

Mode  Identity 

TFF-UM 

TFF-RM 

EFF 

Expt. 

(Approx.) 

44 

42.2 

128.3 

128.7 

Tz 

112 

109.7 

106.6 

107.6 

Stretch 

21 

36.4 

36.5 

36.4 

Bendl 

125 

135.2 

31.4 

31.9 

Txy  1 

15 

21.6 

17.7 

17.7 

17. 9b 

Bend2 

48 

46.0 

6.7 

7.6 

X) 

oo 
o d 

Txy2 

(?) 

(6.95) 

(0.06) 

(0.06) 

(Rot.) 

(?) 

(0.46) 

(0.00) 

(0.00) 

(Tr.) 

C6H6(CH4)i 


Ref.  27 

Eq.  1 

Eq.  1 

Eq.  1 

Mode  Identity 

TFF-UM 

TFF-RM 

EFF 

Expt. 

(Approx.) 

89 

93.4 

90.9 

91.1 

Tx 

89 

93.4 

90.9 

91.1 

Ty 

82 

82.7 

82.7 

82.8 

Stretch 

1  6 

26.8 

25.9 

26.0 

27. 3a 

Bend  y 

16 

26.8 

25.9 

26.0 

27. 3a 

Bend  x 

28 

32.1 

11.8 

11.9 

16. la 

Tz 

(?) 

(6.44) 

(0.63) 

(0.02) 

(Rot.) 

(?) 

(0.41) 

(0.00) 

(0.00) 

(Tr.) 

aIn  the  6q  region. 
bIn  the  origin  region. 


Figure  Captions 


Figure  1.  Water  dimer  geometry,  a.  Geometry  of  the  minimum  energy 
configuration  obtained  from  Scheraga's  intermolecular  potential;16  b. 
Geometry  of  the  minimum  energy  configuration  obtained  from  Jorgensen's 
intermolecular  potential;21  c.  Geometry  at  an  energy  saddle  point  obtained 
from  Scheraga's  intermolecular  potential.16 

Figure  2.  A  few  geometries  of  minimum  energy  configuration  of  higher 
order  water  clusters  obtained  from  Scheraga's  intermolecular  potential.16 
Structures  on  the  left  are  generated  by  the  energy  minimization  routine. 
Structures  on  the  right  are  generated  by  the  TFF-RM  algorithm.  a.  one  of 

the  water  trimers;  b.  one  of  the  water  tetramers;  c.  one  of  the  water 
pentamers;  d.  one  of  the  water  hexamers;  e.  one  of  the  water  heptamers. 

Figure  3.  Cluster  geometries  of  minimum  energy  configuration,  a. 
Benzene( water)]  from  Scheraga's  potential;16  b.  One  of  the  benzene  (water)i 
structures  employing  Jonsson's  potential;23  c.  One  of  the  benzene(ammonia)] 
structures  employing  Scheraga's  potential;16  d.  Second  benzene(ammonia)] 
employing  Scheraga’s  potential;16  e.  Benzene(methane)]  employing 
Scheraga's  potential;16  f.  Second  benzene( water)]  structure  employing 
Jonsson's  potential.23 

Figure  4.  Normal  vdW  mode  vibrations  of  the  benzene( water)]  cluster,  a. 
TFF-RM  method  employing  Scheraga's  potential;16  b.  EFF  method  employing 
Scheraga's  potential;16  c.  EFF  method  employing  Jonsson's  potential.23  The 
mode  energies  are  indicated  for  each  pattern  displayed  (in  cm'1). 
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